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Observations of eclipses of the radio pulsar B1259-63 by the disk of its Be-star companion SS 2883 provide an 
excellent opportunity to study the winds of stars of this type. The eclipses lead to variations in the radio flux 
(due to variations in the free-free absorption), dispersion measure, rotation measure, and linear polarization 
of the pulsar. We have carried out numerical modeling of the parameters of the Be-star wind and compared 
the results with observations. The analysis assumes that the Be-star wind has two components: a disk wind 
in the equatorial plane of the Be star with a power-law fall-off in the electron density n e with distance from 
the center of the star p (n e ~ p _/3 °), and a spherical wind above the poles. The parameters for a disk model 
of the wind are estimated. The disk is thin (opening angle 7.5°) and dense (electron density at the stellar 
surface no e ~ 10 12 /cm 3 , /3 = 2.55), The spherical wind is weak (no e < 10 9 /cm 3 , j3o = 2). 



1 Introduction 

The first binary system consisting of a radio pul- 
sar and a Be star was discovered more than ten 
years ago using the Parkes radio telescope [TJ [2] . 
This is the B1259-63 system, in which the pul- 
sar moves along a very elongated orbit (e > 0.87) 
around its companion - the 10 m Be star SS 2883. 

The first theoretical estimates of the number of 
such systems, carried out in 1983-1987 using the 
"Scenario Machine" [3]-[S] showed that approxi- 
mately one in 700 observed radio pulsars should 
have an OB companion. The detection of the 
B1259-63 system not only confirmed the possiblity 
of the evolutionary scenario obtained by the au- 
thors of [3]-[S], but also provided a powerful tool for 
studies of the associated stellar wind, as was also 
predicted in these studies. Currently, 1500 radio 
pulsars are known, of which at least two are in bi- 
nary systems with OB companions (B1259-63 and 
J0045-7319), consistent with the above estimate. 

In 1997-1998, possible evolutionary tracks for the 
B1259-63 and J0045-7319 systems were computed 
using the "Scenario Machine" [?]• These computa- 
tions were based on an evolutionary scenario that 
predicts the existence of systems consisting of a ra- 
dio pulsar and a massive optical component. The 
high eccentricity of the orbit is explained as an ef- 
fect of the "kick" given by the anisotropic super- 
nova explosion that gave rise to the pulsar. Possible 
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magnitudes and the direction of the kick velocity in 
the B1259-63 system were also derived [5]. 

In 2003, the "Scenario Machine" was used to es- 
timate the number of such systems in which the 
orbital plane of the pulsar and the equatorial plane 
of the Be star do not coincide For characteris- 
tic anisotropic kick velocities of 50-200 km/s, these 
two planes should be inclined relative to one an- 
other by several tens of degrees in more than half 
of these systems. 

A Be star is a main-sequence star of spectral class 
B which has one or more Balmer emission lines 
in its spectrum, with these lines usually display- 
ing two peaks. Struve proposed in 1931 that this 
spectral characteristic could be explained as radi- 
ation from a rotating disk associated with the Be 
star. These disks have now been observed in the 
optical, infrared, and radio. They are comprised of 
dense, slowly rotating material that is located in or 
near the equatorial plane of the Be star. In addi- 
tion, Be stars produce winds with low density and 
high velocity. 

Observations of the transit of PSR B1259-63 pro- 
vide a unique opportunity for studying the charac- 
teristics of the disk of the Be star based on the 
observed variations of the radio flux, linear polar- 
ization, rotation measure, and pulse delay. This is 
possible because the disk of SS 2883 is inclined rel- 
ative to the orbital plane of its companion, so that 
the pulsar is sometimes eclipsed by the disk. 

Attempts were undertaken to construct a model 
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for the Be-star wind and compare the calculated 
parameters with the observed dispersion measure, 
rotation measure, and pulse delay time. A disk 
model with an exponential fall-off of the electron 
density with distance from the Be star and with 
height above the plane of the disk was considered 
in [IU]. 

A model with a power-law fall-off of the den- 
sity in the disk is more physically justified, how- 
ever [111; with some parameters, this model was 
able to explain the observed variations in the dis- 
persion measure near periastron. A model with a 
clumpy disk wind from the Be star has also been 
developed [12 [T3] (in addition to, not replacing 
the power-law disk). Rapidly moving bubbles (ve- 
locities of « 2000 km/s) with electron densities 
that differ from the value in the surrounding re- 
gion (n e ~ 10 6 /cm 3 , and dimensions < 10 10 cm at 
a distance of 20-50 stellar radii) give rise to appre- 
ciable fluctuations in the electron density along the 
line of sight, and thereby to fluctuations in the flux, 
dispersion measure, and rotation measure. 

An attempt to establish the position of the or- 
bital plane of PSR B1259-63 relative to the equa- 
torial plane of SS 2883 was made using timing 
measurements for the pulsar [Ml [15]. However, it 
proved impossible to construct a unique model for 
the system due to the noise level and the prolonged 
eclipses of the pulsating radio emission. 

A disk model for the Be-star wind is also con- 
sidered, for example, in [Trj],ll7|. where the ratio N 
of the densities of the outflowing material in the 
equatorial plane and at the poles is derived. If 
the Be star rotates with a velocity v that is 90% 
of the critical velocity v cr a then N = 150; slower 
rotation lowers this value (when v = 70% • v cr u, 
N w 15 — 20). According to [18], the rotational 
velocity of SS 2883 is approximately 70% of the 
critical velocity. 

We consider here the absorption of the radio 
emission in a model with a power-law disk wind 
in the equatorial plane and a spherical wind above 
the poles of the Be star (without including the ef- 
fect of dumpiness). In contrast to previous studies, 
we compare the calculated and observed fluxes of 
pulsating radio emission. The derived parameters 
of the wind differ from those presented in [IT] . The 
available rotation measure data are sufficient only 
to derive order of magnitude estimates of the mag- 
netic field strength in the wind (in particular, in 




Fig. 1 : Schematic model for a system consisting of 
a radio pulsar and a Be star. The X axis is directed 
along the line of sight away form the observer, uj is 
the longitude of the ascending node, i is the orbital 
inclination (the angle between the plane of the sky 
YZ and the orbital plane of the pulsar), P indicates 
the pulsar and Be the Be star, it, is the angle be- 
tween the XY plane and the plane of the Be-star 
disk, and Of, is the angle between the Y axis and 
the line corresponding to the intersection of the XY 
plane and the disk of the Be star. 

bubbles), as was done in [UJ-[I3]. Accordingly, we 
will consider further only the radio flux and disper- 
sion measure. 

2 Model Be-star wind 

Figure [JJ depicts a schematic model of a binary 
system consisting of a radio pulsar and a Be star. 
In this figure, u> is the longitude of the ascending 
node, i is the orbital inclination (the angle between 
the plane of the sky YZ and the orbital plane of 
the pulsar), P indicates the pulsar and Be the Be 
star, if, is the angle between the XY plane and the 
plane of the Be-star disk, and fit, is the angle be- 
tween the Y axis and the line corresponding to the 
intersection of the XY plane and the disk of the Be 
star. The X axis is oriented along the line of sight 
away from the observer. 

The calculated model for the wind of SS 2883 is 
presented in Fig. [5] There is a disk wind with open- 
ing angle ip in the equatorial plane, and a spheri- 
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Fig. 2: Model Be-star wind. There is a disk wind 
with a power-law fall-off in the electron density in 
the equatorial plane and a spherical wind above the 
poles. 



cal wind above the poles, which we take here to be 
isothermal and to have a constant radial speed. We 
assume that the outflowing gas is fully ionized. 

The electron density n e a distance p from the 
star is calculated using the formula 



n 0e — 
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Here, no e is the electron density at the stellar sur- 
face (i.e., at a distance po from its center). For the 
spherical wind, @q — 2; this is a free parameter for 
the disk wind, whose value is one of the parameters 
for which we are searching. 

We also assumed that the temperature of the disk 
wind can be determined by the power-law relation 
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Here, T is the temperature of the plasma at a dis- 
tance p from the star and To is the temperature at 
the stellar surface. 

The optical depth to radio emission can be cal- 
culated using the formula [19] : 
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Here, v is the frequency of the radio emission, T the 
temperature of the wind (J5J), and EM the emission 
measure, which is defined to be 

EM = J nidi. (4) 
The dispersion measure is given by the formula 

DM = I n R dl. (5) 
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Fig. 3: Dependence of the 1384 MHz radio flux of 
B1259-63 on the position of the pulsar relative to 
SS 2883; v is the true anomaly. The hollow squares 
show the observational data with their errors. The 
solid curve denotes the calculated model. 

The quantity n e in ((4| and ([5|) is calculated using 
©■ 



3 Results 

We carried out theoretical calculations of the flux 
and dispersion measure of the radio emission of the 
B1259-63 pulsar as a function of its position rela- 
tive to its companion - the Be star SS 2883. These 
calculated values were compared with the observa- 
tions. 

The observational data on the radio fluxes were 
taken from Table 2 of [13], and on the dispersion 
measures from Table 1 of [TU], Table 2 of [TT], and 
Table 3 of [T3]. The longitude of the ascending 
node of the pulsar orbit in the B1259-63 system 
is oj — 138° PJ, the orbital inclination is i = 36°, 
the characteristic radius of the star is po = 6i?© 
the mass of the optical star is M = 10 M Q and 
the projection of the semi-major axis onto the line 
of sight is asini= 1295.98 light seconds [2J. The 
interstellar contribution to the dispersion measure 
of the B1259-63 pulsar is 146.8 Pc/cm 3 [15] . 

The velocity of the Be-star wind at infinity is 
V^o 1000 — 1300 km/s in the equatorial plane and, 
w 2300 km/s at the poles [T7] (the characteristic 
wind velocity at infinity for SS 2883 is « 1350 ±200 
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Fig. 4: Same as Fig. 3 for 2496 MHz. 

km/s [2H]). The radial outflow velocity of the disk 
material at the Be-star surface is 5-10 km/s |21j . 

To check the resulting parameters of the stellar 
wind, we must compare the mass-loss rate in the 
model obtained with the characteristic values for 
Be stars. The rate at which matter flows from the 
star can be found from the expression 

M = Clp 2 nmiV; (6) 

where p is the distance from the star, n the parti- 
cle density of the wind at the distance p which is 
calculated using (fTJ), rtii the mean mass of the ions 
in the wind, V the wind velocity at the distance p, 
and Q the solid angle into which the given type of 
wind (spherical or disk) flows. 

The mass loss of Be stars was investigated in de- 
tail in 22 . A Be star with an initial mass of 10M Q 
loses mass at the rate M = 3 • 10~ 9 M Q /year. 

We carried out numerical computations in order 
to determine the best-fit orientation of the disk of 
the Be star SS 2883 relative to the orbital plane of 
PSR B1259-63, which is given by the parameters fib 
and if, (Fig. [T]). We also wished to determine the 
best-fit values of fto for the disk wind and the elec- 
tron density at the stellar surface no for both the 
disk and spherical winds (/3q = 2 for the spherical 
wind). We constructed theoretical light curves for 
the pulsar at 1384, 2496, 4800, and 8400 MHz to- 
gether with the dependence of the dispersion mea- 
sure of the pulsar radio emission on the position of 
the pulsar relative to its companion. 



Fig. 5: Same as Fig. 3 for 4800 MHz. 

The results of the numerical simulations are 
shown in[3]-[7]together with the observational data. 
The criterion used to evaluate the correctness of a 
model is the closeness of the computed curves to the 
observational data. In addition, the model should 
not contradict our understanding of the matter out- 
flow rate in Be stars. 

The computations yielded the following best-fit 
parameters for the disk model. The disk opening 
angle was tp — 7.5°, the electron density in the 
disk at the stellar surface was no e w 1 • 10 12 /cm 3 , 
and fla = 2.55. The spherical wind was weak 
(«0e < 10 9 /cm 3 , fa = 2). The orientation of the SS 
2883 disk relative to the orbital plane is described 
by the parameters and it, (Fig. [T]), which have 
the values fib = 12° and % — 67° in the best-fit 
model. The disappearance of the pulsating emis- 
sion of B1259-63 near periastron is due to eclipses 
of the pulsar by the Be-star disk. The weak spheri- 
cal wind does not make an appreciable contribution 
to either the decrease in the radio flux or the in- 
crease in the dispersion measure. The model mass- 
loss rate via the SS 2883 wind estimated using Q 
is consistent with our current understanding of the 
mass-loss rates of Be stars: M w 3- 10~ 9 Af Q /yr, in 
good agreement with the results of [2"2"] . 

Although we have assumed that the temperature 
of the wind may depend on the distance from the 
Be star ([2]), the best- fit model was obtained for a 
constant temperature, equal to 10 4 K. Substituting 
([T|) and into we can see that the behavior of 
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Fig. 6: Same as Fig. 3 for 8400 MHz. 



the light curve will not change if 2(3q — 1.35/3i = 5.1 
if we suppose that (5\ / (so that consequently, 
(3q 7^ 2.55). However, the variations of (3$ must not 
be in contradiction with the observed dispersion 
measures, and cannot be large. 



Fig. 7: Dependence of the dispersion measure on 
the true anomaly v of the B1259-63 pulsar. The 
hollow squares indicate the observational data with 
their errors. The solid curve denotes the model 
with a thin and dense disk wind and a negligi- 
bly small spherical wind (we took the interstellar 
contribution to the dispersion measure to be 146.8 
pc/cm3). 



Figure [3] clearly displays a secondary minimum 
of the light curve at values of the true anomaly v 
from 100° to 120°. Our model for the wind is not 
able to fit this observational feature. The compu- 
tations yielded a model with flf, w —55° that can 
fit this secondary minimum, whose presence is then 
explained by the absorption of the radio emission 
in the disk. However, in this case, the pulsar and 
observer are located on the same side of the disk at 
periastron, and we must introduce a strong spher- 
ical wind in order to fit the decrease in the radio 
flux. This led to an implausibly high mass-loss rate 
(M ss 3- 10 _7 M©/yr). In addition, a strong spheri- 
cal wind would make an appreciable contribution to 
the dispersion measure, so that it would be neces- 
sary to appreciably lower the interstellar dispersion 
measure (from 146.8 PC/cm 3 to 135-140 Pc/cm 3 ). 
We, accordingly, rejected this second model. 

The parameters of the best-fit model presented 
here differ substantially from the results of [TT| . 
where, in particular, it was found that (3q — 4.2 in 
the disk and that the mass-loss rate of the optical 
star is M = 5 ■ 10~ 8 Af o /yr. 



4 Conclusion 



Our computations have yielded the following best 
fit paramters for the disk model used. The opening 
angle of the disk is <p = 7.5°, the electron density in 
the disk at the stellar surface is n^ e w 1 • 10 12 /cm 3 , 
(3q = 2.55. The spherical wind is weak (no e < 
10 9 /cm 3 , /3o = 2) and cannot make a significant 
contribution to either the decrease in the radio flux 
or the increase in the dispersion measure near peri- 
astron. The orientation of the SS 2883 disk relative 
to the pulsars orbital plane is fit by the parameters 
n b = 12° and i b = 67° (Fig. [[]). The mass-loss 
rate in the SS 2883 wind estimated using ©, is 
consistent with current thinking about the mass- 
loss rates of Be stars: M « 3 • 10 _9 MQ/yr, in good 
agreement with the results of [52]. The disappear- 
ance of the pulsating radio emission of B1259-63 
near periastron is due to eclipsing of the pulsar by 
the dense disk of the Be star. 

The author thanks the anonymous referee for 
useful comments. 
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